Relativistic many-body calculations of excitation energies and transition rates in 

ytterbium-like ions 
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Excitation energies, oscillator strengths, and transition rates are calculated for {5cP + 5d6s + 6s^)- 
{5d6p + 5d5/ + 6s6p) electric dipole transitions in Yb-like ions with nuclear charges Z ranging from 
72 to 100. Relativistic many-body perturbation theory (RMBPT), including the retarded Breit 
interaction, is used to evaluate retarded El matrix elements in length and velocity forms. The 
calculations start from a [Xe]4/^* core Dirac-Fock potential. First-order RMBPT is used to obtain 
intermediate coupling coefficients, and second-order RMBPT is used to determine matrix elements. 
A detailed discussion of the various contributions to energy levels and dipole matrix elements is given 
for ytterbium like rhenium, Z=75. The resulting transition energies are compared with experimental 
values and with results from other recent calculations. Trends of excitation energies, line strengths, 
oscillator strengths, and transition rates as functions of nuclear charge Z axe shown graphically for 
selected states and transitions. These calculations are presented as a theoretical benchmark for 
comparison with experiment and theory. 

PACS numbers: PACS: 32.70.Cs, 31.15.Md, 31.25.Eb, 31.25.Jf, 31.30.Jv 



I. INTRODUCTION 

We report results of ab initio calculations of excitation 
energies, oscillator strengths, and transition rates in Yb- 
like ions with nuclear charges Z ranging from 72 to 100. 
The ions considered here, starting from doubly ionized 
Hf III, all have bcP ground states. We do not consider 
Yb I and Lu II both of which have a 6s^ ground state con- 
figuration. In recent publications |l], || ||, ^, ^, |^, |^ , the 
spectra of Re VI, Os VII, and Ir VIII were studied and 
energies levels of the 5cP, 5d6s, 5d6p, 5d5f and 6s6p con- 
figurations were determined. The Cowan Hartree-Fock 
code 1^ with relativistic and correlation options was used 
in Refs. 0, ||, ||, |4[ |[ |[ Q to calculate energy levels and 
to carry out least-squares adjustments of energies. 

Although we do not consider Yb I and Lu II here, it 
should be noted that Porsev et al. Q recently carried 
out elaborate calculations of electric-dipole amplitudes 
in atomic ytterbium. Moreover, Martin et al. |10 listed 
energies for 249 levels of Yb I and 40 levels of Lu II. 

In the present paper, we use relativistic many-body 
perturbation theory (RMBPT) to determine energies of 
the 14 even-parity 5d^, 5d6s, and 6s^ states and the 
36 odd-parity 5d6p, 5d5f, and 6s6p states for Yb-like 
ions. We illustrate our calculation with a detailed study 
of Re VI, Z= 75. Our first-order RMBPT calculations 



include both the Coulomb and retarded Breit interac- 
tions but our second-order calculations are limited to the 
Coulomb interaction only. 

Reduced matrix elements, line strengths, oscillator 
strengths, and transition rates are determined for all al- 
lowed and forbidden electric-dipole transitions between 
even-parity {5d^ + 5d6s + 6s^) and odd-parity {5d6p + 
5d5f + 6s6p) states. Retarded El matrix elements are 
evaluated in both length and velocity forms. RMBPT 
calculations that start from a local potential are gauge 
independent order-by-order, provided "derivative terms" 
are included in second- and higher-order matrix elements 
and careful attention is paid to negative-energy states. 
The present calculations start from a nonlocal [Xe]4/^^ 
Dirac-Fock (DF) potential and consequently give gauge- 
dependent transition matrix elements. Second-order cor- 
relation corrections compensate almost exactly for the 
gauge dependence of the first-order matrix elements and 
lead to corrected matrix elements that differ by less than 
5% in length and velocity forms for all of the ions con- 
sidered here. 

Energies from the present calculation agree well with 
results given in Refs. ||l|, ^ ||, |[ ||, ||, 0] for low-lying lev- 
els, but disagree substantially for various highly-excited 
levels, as discussed later. 
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The RMBPT formalism developed previously ||Tl|, |T|, 
[l^ , |l^ , [l5|| for Be-, Mg-, and Ca-like ions is used here to 
describe perturbed wave functions, to obtain the second- 
order energies Um, and to evaluate first- and second- 
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TABLE I: Possible two-particle states in the Yb-like ions 



jj coupling 


LS coupling 


jj coupling 


LS coupling 


5d3/25d3/2 


(0) 




Po 


5d3/26pi/2(2) 


5d6p 


'F2 




(0) 




So 


5d3/26p3/2(2) 


5d6p 


'D2 


6S1/26S1/2 


0) 




So 


5d5/26pi/2(2) 


5d6p 


'D2 










5d5/26p3/2(2) 


5d6p 


'P2 


5^3/25^5/2 


(1) 


5d^ ^ 


Pi 


5d3/25/5/2(2) 


6s6p 


'P2 


5d3/26Si/2 


;i) 


5d6s 




5d3/25/7/2(2) 


5d5/ 


'F2 










5d5/25/5/2(2) 


5d5/ 


'D2 


5d3/25d3/2 


(2) 




F2 


5d5/25/7/2(2) 


5d5/ 


'D2 


5^3/25^5/2 


(2) 




D2 


6Sl/26p3/2(2) 


5d5/ 


'P2 


5d5/25d5/2 


(2) 




P2 








5d3/26Si/2 


;2) 


5d6s 




5d3/26p3/2(3) 


5d6p 


'Fs 


5d5/26Si/2 


;2) 


5d6s 




5d5/26pi/2(3) 


5d6p 


'D, 










5d5/26p3/2(3) 


5d6p 


'Fs 


5d3/25d5/2 


(3) 


5d^ ^ 


F3 


5d3/25/5/2(3) 


5d5/ 


'Fs 


5d5/26Si/2 


;3) 


5d6s 




5d3/25/7/2(3) 


5d5/ 


'G^ 










5d5/25/5/2(3) 


5d5/ 


'Ds 


5^3/25^5/2 


(4) 


5d^ ^ 


Fa 


5d5/25/7/2(3) 


5d5/ 


'Fs 


5d5/25d5/2 


(4) 




















5d5/2 6p3/2(4) 


5d6p 




5^3/26^3/2 


(0) 


5d6p 


'Po 


5d3/25/5/2(4) 


5d5/ 




5^5/25/5/2 


;o) 


6s6p 


'Po 


5d3/2 5/7/2(4) 


5d5/ 


'H, 


6S1/26P1/2 


0) 


5d5/ 


'Po 


5d5/25/5/2(4) 


5d5/ 


'F, 










5d5/2 5/7/2(4) 


5d5/ 




5^3/26^1/2 


(1) 


5d6p 


3 7-» 








5^3/26^3/2 


(1) 


5d6p 


'Pi 


5d3/9 5f7/2(5) 


5d5/ 


'H, 


5^5/26^3/2 


(1) 


5d6p 


'Pi 


5d5/25/5/2(5) 


5d5/ 


'G, 


5^3/25/5/2 


Jl) 


6s6p 


'Pi 


5d5/25/7/2(5) 


5d5/ 


'H, 


5^5/25/5/2 


[1) 


5d5/ 


'D, 








5^5/25/7/2 


Jl) 


6s6p 


'Pi 


5d5/25/7/2(6) 


5d5/ 


'H, 


6S1/26P1/2 


1) 


5d5/ 


'Pi 








6S1/26P3/2 


1) 


5d5/ 


'Pi 









TABLE IL Contributions to energy matrices E[Q, Q'] 

(a.u.) for odd-parity states Q—nilijin2l2j2{J), Q'= 
n^hj^n^lAjA^J) with J=Q before diagonalization in the case 
of Yb-like rhenium, Z — 75. 

Q, Q' g(^) B(^) 

5d3/26p3/2, 5d3/26p3/2 -3.99398 0.42731 0.00679 -0.14018 

5d5/25/5/2, 5d5/25/5/2 -3.42978 0.44889 0.00430 -0.16664 

6si/26pi/2, 6si/26pi/2 -3.64181 0.29737 0.00591 -0.10004 

5d3/26p3/2, 5d5/25/5/2 -0.01871 0.00003 0.01465 

5d5/25/5/2, 5d3/26p3/2 -0.01871 0.00003 0.00965 

5d3/26p3/2, 6S1/26P1/2 0.02191 -0.00001 -0.01278 

6S1/26P1/2, 5d3/26p3/2 0.02191 -0.00001 -0.00996 

5d5/25/5/2, 6S1/26P1/2 -0.01992 -0.00001 0.00221 

6S1/26P1/2, 5d5/25/5/2 -0.01992 -0.00001 0.00252 



order transition matrix elements |13[| . Ions of the Yb 
sequence, starting from the Hf III ion | [T^ and continu- 
ing onward have a 5(P ground state. This is similar to 
the previously studied Ca sequence ||l^ , where ions start- 
ing from Ti III have a 3(fi ground state. The primary 
differences between calculations for Ca-like and Yb-like 
ions arise from the increased number of the orbitals in 
the DF core potential, [Xe]4/i'* instead of [Ar] ([Ar] = 



Is22s22p63s23p6 and [Xe] = [AT]3d^°4s^4p^Ad''°5s^5p'^), 
and the strong mixing between both even-parity {5(P + 
5d6s + Qs"^) and odd-parity {5d6p + 5d5f + 6s6p) states. 
These differences lead to much more laborious numeri- 
cal calculations. The calculations are carried out using 
sets of DF basis orbitals that are linear combinations of 
B-splines. These B-spline basis orbitals are determined 
using the method described in Ref. We use 40 B- 

splines of order 8 for each single-particle angular momen- 
tum state and we include all orbitals with orbital angular 
momentum Z < 7 in our basis set. 



A. Model space 

The model spaces for the (Sd^ -|- 5d6s + 6s^) and 
{5d6p + 5d5/ -I- 6s6p) complexes in Yb-like ions have 
14 even-parity states and 36 odd-parity states, respec- 
tively. These states are summarized in Table |, where 
both jj and LS designations are given. When start- 
ing calculations from DF wave functions, it is natural 
to use jj designations for uncoupled matrix elements; 
however, neither jj- nor iS'-coupling describes physical 
states properly, except for the single-configuration state 
5^5/25/7/2 (6) = 5d5f '^Hq. The strong mixing between 
5d6p, 5d5f, and 6s6p states was discussed previously in 
Refs. 

B. Example: energy matrix for Re^® 

Details of the theoretical method used to evaluate 
second-order energies for ions with two valence electrons 
are given in Refs. (l^, |l^ and will not be repeated 
here. The energy calculations are illustrated in Table 
U, where we list contributions to the energies of odd- 
parity J = states of Re+^. We present zeroth-, first-, 
and second-order Coulomb energies 
together with the first-order retarded Breit corrections 
B^^^ [|l8j. It should be mentioned that the difference be- 
tween first-order Breit corrections calculated with and 
without retardation is less than 2%. As one can see 
from Table |l[ the ratio of off-diagonal to diagonal ma- 
trix elements is larger for second-order contributions than 
for first-order contributions. Another difference between 
first- and second-order contributions concerns symme- 
try properties: first-order off-diagonal matrix elements 
are symmetric, whereas second-order off-diagonal matrix 
elements are unsymmetric (Lindgren and Morrison |l9| , 
chap. 9]). Indeed, E'^^'>[Q,Q'] and E'^^^[Q',Q] differ in 
some cases by 20-50% and occasionally even have oppo- 
site signs. The ratio of off-diagonal to diagonal matrix 
elements for Breit corrections i?*-^) is much smaller than 
for Coulomb corrections. 

After evaluating the energy matrices, we calculate 
eigenvalues and eigenvectors for states with given val- 
ues of J and parity. There are two possible methods 
to carry out the diagonalization: either diagonalize the 
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TABLE III: Energy of 5<r , 5d6s, 6s\ 
and 6s6p states in Yb-like Re"""^ (cm~') 



5d6p, 5d5/, 
Notation: 





£■(0+1) 




^(tot) 


^(0+1) 




^(cxc) 


Level 


Absolute energies 


Excitation energies 




-1184024 


-40817 


-1224841 











5d^ 


-1177178 


-39172 


-1216350 


6847 


1645 


8491 


5d^ ^Fi 


-1170966 


-38780 


-1209746 


13058 


2037 


15095 




-1167037 


-43309 


-1210346 


16987 


-2493 


14495 


5d^ 


-1166402 


-41572 


-1207975 


17622 


-756 


16866 


5d^ -^Pi 


-1163473 


-41812 


-1205285 


20551 


-995 


19556 


5d^ ^P2 


-1156318 


-40152 


-1196470 


27706 


665 


28371 


5d^ 


-1155966 


-42068 


-1198034 


28059 


-1251 


26807 


5d^ 'So 


-1127568 


-46217 


-1173785 


56456 


-5400 


51056 


5d6s ^Di 


-1097044 


-35141 


-1132185 


86980 


5675 


92656 


5d6s ^D2 


-1094845 


-35442 


-1130287 


89180 


5375 


94554 


5d6s ^Ds 


-1087729 


-33669 


-1121397 


96296 


7148 


103444 


5d6s 'l>2 


-1079142 


-36325 


-1115467 


104883 


4491 


109374 


5d6p 


-1030357 


-30866 


-1061223 


153667 


9951 


163618 


5d6p ^L>i 


-1024847 


-33457 


-1058303 


159178 


7360 


166538 


5d6p ^D2 


-1018155 


-30518 


-1048673 


165869 


10299 


176168 


5d6p ^Fs 


-1017295 


-30334 


-1047629 


166729 


10483 


177212 


5d6p 


-1012677 


-29311 


-1041988 


171347 


11505 


182853 


5d6p 


-1007838 


-30332 


-1038169 


176187 


10485 


186672 


5d6p ^Pi 


-1005870 


-31934 


-1037804 


178154 


8882 


187037 


5d6p ^Po 


-1004380 


-31262 


-1035642 


179644 


9554 


189199 


5d6p 


-1001386 


-28009 


-1029395 


182639 


12807 


195446 


5d6p ^P2 


-999246 


-29109 


-1028355 


184778 


11708 


196486 


5d6p 


-995855 


-31244 


-1027099 


188169 


9573 


197742 




-991963 


-34225 


-1026189 


192061 


6592 


198652 


5d6p ^Pi 


-989455 


-34880 


-1024335 


194569 


5936 


200506 


6s6p ^Po 


-941422 


-21782 


-963204 


242603 


19035 


261637 


6s6p ^Pi 


-936597 


-22693 


-959291 


247427 


18123 


265550 


6s6p ^P2 


-921614 


-21582 


-943197 


262410 


19234 


281644 


5d5/ 'G4 


-914890 


-27340 


-942229 


269135 


13477 


282611 


5d5/ ='//4 


-911078 


-28960 


-940038 


272946 


11857 


284803 


5d5/ ^//s 


-910564 


-26917 


-937481 


273460 


13900 


287360 


5d5/ 3P2 


-909842 


-32013 


-941855 


274182 


8804 


282986 


5d5/ ^Pa 


-907828 


-33638 


-941466 


276197 


7178 


283375 


5d5/ ^'//e 


-905206 


-23491 


-928698 


278818 


17325 


296143 


5d5/ ^P4 


-903304 


-28769 


-932074 


280720 


12047 


292767 


5d5/ '732 


-900678 


-34598 


-935275 


283347 


6219 


289565 


5d5/ ^Di 


-898884 


-36856 


-935740 


285140 


3960 


289100 


5d5/ ^^Gs 


-897436 


-41168 


-938604 


286588 


-352 


286237 


5d5/ =^P»2 


-893563 


-36014 


-929576 


290462 


4803 


295265 


6s6p 'Pi 


-892813 


-36571 


-929383 


291211 


4246 


295457 


5d5/ ^G4 


-892131 


-38708 


-930839 


291893 


2108 


294002 


5d5/ ^P»3 


-891083 


-38081 


-929165 


292941 


2735 


295676 


5d5/ ^Gs 


-889833 


-39291 


-929124 


294191 


1526 


295717 


5d5/ ='P2 


-887212 


-36551 


-923763 


296812 


4266 


301078 


5d5/ =^Pi 


-886406 


-36865 


-923271 


297619 


3952 


301570 


5d5/ ^Po 


-886124 


-36385 


-922510 


297900 


4431 


302331 


5d5/ 'P3 


-884984 


-39786 


-924770 


299040 


1031 


300071 


5d5/ '1/5 


-872058 


-49224 


-921282 


311966 


-8407 


303559 


5d5/ 'Pi 


-866971 


-40448 


-907419 


317054 


369 


317422 



sum of zeroth- and first-order matrices, then calculate 
the second-order contributions using the resulting eigen- 
vectors; or diagonalize the sum of the zeroth-, first- and 
second-order matrices together. Following Ref. we 
choose the second method here. 



The energy calculations are illustrated for Re+^ in Ta- 
ble III. Energies listed under the heading "Absolute 
energies" are given relative to the [Xe]4/^^ core, while 
those listed under the heading "Excitation energies" are 
given relative to the 5c?^ ground state. In the ta- 



ble, we present absolute zeroth- plus first-order Coulomb 
and Breit energies £;("+') = -|- -f S^), absolute 
second-order Coulomb energies 

£;(2) ^ and the sum 

^(tot)^ 

We also give the breakdown excitation energies and total 
excitation energies E^'^^''\ As can be seen from the table, 
the second-order contribution is about 3% of the abso- 
lute energy but accounts for 5% - 15% of the excitation 
energy. This table clearly illustrates the importance of 
including second-order contributions. As mentioned pre- 
viously, neither jj- nor LS'-coupling describes physical 
states properly; nevertheless, we use LS designations to 
label levels. We organize levels in the table according to 
decreasing values of E^^~^^\ After including second-order 
corrections, the present ordering differs in some cases 
from the ordering according to decreasing 

^(tot)^ Thus, 

the ordering of 5cP and Sd^ ^Pq levels, 5d5/ ^H^ and 
bdbf levels, and 5d5/ and 6s6p ^Pi levels are 
interchanged in Re"*"^ after second-order corrections are 
added. 



Problems arising when using different model spaces in 
RMBPT theory were examined by Johnson et al. ||l2|. A 
major difference between Yb-like and Be-like systems, is 
that we could not construct as complete a model space 
for a two-electron system with a|Xe]4/^'* core as we did 
for a system with a [He] core To do so would re- 

quire us to include all possible two-particle states that 
could be constructed from unoccupied n=5 and n—6 or- 
bitals in our model space. In the case of Be-like ions, 
the model space is much simpler, being constructed from 
n=2 orbitals only. A second, but related, problem is that 
different model spaces lead to different results. For ex- 
ample, energy levels of Hg I, calculated using RMBPT 
with (6s6p) and (6s^ -I- 6p^) model spaces differed from 
those calculated with {6s6p+6p6d) and (6s^-f 6p^-|-6s6(i) 
model spaces by about 500 cm^^ jl^. We confirm this 
result here comparing calculations of 5p6d energy levels 
starting from a 5p6d model space and those starting from 
{5p6d + 5d5/) or (5p6d + 5rf5/ + 6s6p) model spaces in 
Re"*"^. The largest difference occurs in the E^'^\ which 
changes by 1000 cm~^ in some cases. The change in 
^(o-(-i) jg snialler by a factor of two and has an opposite 
sign. The resulting change in 75^*°*^ is about 500 cm~^. 
Similar tests of model-space dependence along the Yb se- 
quence set a limit of about 500 cm~^ on the accuracy of 
the present second-order calculations. 
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C. Z-dependence of energies 



compensated by "derivative terms" p(dcrv) 



One unique feature of the present calculations is 
the inclusion of second-order correlation corrections. 
We illustrate the Z-dependence of the second-order 
energy E^^'' in Fig. |] for even-parity levels with 
J=2 {5(f 3^2, ^D2, and 5d6s ^'^Da) and J=3, 4 
(5^2 ^Fa, ^Fi, ^Gi and 5d6s ^£>3). As can see from 
this figure, the second-order energy, 

£;(2) slowly increases 
with Z in the range 3-5x10^ cm~^. The smooth Z- 
dependcnce for these nine terms is exceptional and is not 
found for other terms discussed below. 

Excitation energies even-parity and odd- 

parity states relative to the bcP ground state, divided 
by {Z — 65)^, are shown in Figs. || and |[ Both designa- 
tions are shown in these figures: LS for low Z and jj for 
high Z. The variation of the b(P levels with Z shown 
in Fig. Ija. Strong mixing between ^d^/2^d^/2{J) ^-nd 
^d^/2^d^/2{J) states with J=2 or 4 leads to rapid vari- 
ations with Z in the corresponding Grotrian diagrams. 
Note that the ^Z?2 and levels cross between Z=79 
and 80, ^Po and '^^3 levels cross between Z—80 and 81, 
and "^Pi and levels cross between Z=93 and 94. We 
give excitation energies of the other five even-parity 5d6s 
and 6s^ levels in Fig. Two small sharp features in 
the 6s^ level occur at Z=82 and Z=85. The origin of 
these irregularities is discussed in Appendix A. Energies 
of odd-parity levels with J—0 and 4, including mixing of 
5d6p, bdbf and 6s6p states, are given in Fig. ^. The sharp 
features in the curves describing 5d5f and 6s6p states are 
similar to those mentioned above for even-parity states 
and discussed in Appendix A. Avoided level crossings for 
the odd-parity J = levels are seen Fig. |^ near Z=76 
and Z—SA. Similar avoided crossings can be observed for 
the odd-parity complex with J=4 in Fig. [^b. 



D. Z-dependence of matrix elements for 
electric-dipole transitions 

We designate the first-order dipole matrix element by 
Z^^\ the Coulomb correction to the second-order matrix 
element by Z^"^^ , and the second-order Breit correction by 
5(2). The evaluation of Z'^^\ Z'^^\ and B^^) for Yb-like 
ions follows the pattern of the corresponding calculation 
for beryllium- like ions in Refs. |20|. These matrix ele- 
ments are calculated in both length and velocity gauges. 
Differences between length and velocity forms, are illus- 
trated for the uncoupled 5d3/25d5/2(2)-5c?5/26pi/2(3) ma- 
trix element in Fig. ^. The second-order Breit matrix el- 
ement B^^^ is multiplied by a factor of 50 in order to put 
it on the same scale as the second-order Coulomb matrix 
element Z'^' . The sharp features have the same origin as 
those in the the second-order energy matrix and are dis- 
cussed in Appendix B. Contributions of the second-order 
matrix elements Z^^^ and i?*-^^ are much larger in veloc- 
ity [V) form than in length (L) form as seen in panels 
(a) and (6) of Fig. ^. As shown later, L-V differences are 



E. Example: dipole matrix elements in Re+® 

We list uncoupled first- and second-order dipole ma- 
trix elements Z'^^\ Z'-^\ B'^'^\ together with derivative 
terms p(dcrv) -^^ ^able For simplicity, 



we consider only the 9 dipole transitions between even- 
parity states with J=0 and odd-parity states with J=l. 
The derivative terms shown in Table IV arise because 



transition amplitudes depend on energy, and the transi- 
tion energy changes order-by-order in perturbation the- 
ory. Both L and V forms are given for the matrix el- 
ements. We see that the first-order matrix elements 



Z^l^ and Z\}' differ by 10-20% and that the L-V dif- 
ferences between second-order matrix elements are much 
larger for some transitions. The first-order matrix ele- 
and z{}^ are non-zero for the three matrix 



.(1) 



7W 



5^3/25^3/2 (0)-5d3/26pi/2 (1) , 



5^3/25^3/2(0)- 



elements 

54/26^3/2(1). and 5d5/25d5/2(0)-5(i5/26p3/2(l), but van- 
ish for the remaining six transitions. Moreover, those six 

f2l f21 

transitions, Z)^ and Zy have the same order of mag- 
nitude as the three non-zero first-order matrix elements. 
This confirms the importance of including higher-order 
RMBPT corrections. It can also see from Table 



(dcrv) 
L 
.(1) 



is almost equal to Z 



(1) 



but P, 



(dc 



IV that 



is smaller than 



Zy by five to six orders of magnitude. 

In Table line strengths for Re"^'^ in length and veloc- 
ity forms are given for the J=0-J'=1 transitions consid- 
ered in Table |^as well as for other J- J' transitions. We 
see that L and V forms of the coupled matrix elements in 
Table ^ differ only in the second or third digits. As men- 
tioned in the introduction, these L-V differences arise 
because we start our RMBPT calculations using a non- 
local DF potential. If we were to replace the DF potential 
by a local potential, the differences would disappear com- 
pletely. Another source of L~V differences is the use of 
an incomplete model space. The last two columns in Ta- 
ble show L and V values of line strengths calculated 
without the second-order contribution. As can be seen, 
including second-order corrections significantly decreases 
L-V differences. 



III. RESULTS AND DISCUSSION 

We calculate energies of the 14 even-parity 50?^, 5d6s, 
and Qs^ states as well as the 36 odd-parity 5c?6p, 5c?5/, 
and 6s6p states for Yb-like ions with nuclear charges in 
the range Z = 72 to 100. Reduced matrix elements, line 
strengths, oscillator strengths, and transition rates are 
also determined for all electric-dipole transitions between 
even-parity and odd-parity states for each ion. Compar- 
isons with experimental data and other theoretical results 
are also given. Our results are presented in three parts: 
transition energies, fine-structure energy differences, and 



CO 25 
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- 5d6s 'D^ 
-5d6s'D, 



75 SO 85 90 95 
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-5d^'G, 



75 80 86 90 95 

Nuclear charge Z 



FIG. 1: Z-dependence of the second-order energy E^'^^ for the 5d^ and 5d6s energy levels 
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FIG. 2: Z-dependence of the excitation energy E^*"^"^ /{Z - 65)^ in cm ^ 
for even-parity levels. 
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FIG. 3: Z-dependence of the excitation energy e'-'""'^/{Z - 65)^ in cm ^ 
for odd-parity levels. 
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TABLE IV: Uncoupled reduced matrix elements in length L and velocity V forms for even-odd parity transitions in Rc^ 



even-parity 



odd-parity 



7(1) 



j(l) 

-'V 



7(2) 



7(2) 



B 



(2) 



(2) 



^(dcrv) 



r5(dcrv) 
V 



5d3/25d3/2(0) 5d3/26pi/2(l) 0.84736 0.73795 -0.18176 -0.32201 0.00105 0.00032 0.84733 -0.00001 

5^3/25^3/2(0) 5^3/26^3/2(1) -0.33505 -0.29206 0.00633 0.03976 -0.00074 -0.00053 -0.33502 0.00003 

5d3/25d3/2(0) 5^6/26^3/2(1) 0.00000 0.00000 -0.06416 -0.09997 0.00002 0.00000 0.00000 0.00000 

5d5/25dB/2(0) 5d3/26pi/2(l) 0.00000 0.00000 -0.04075 -0.08423 -0.00001 -0.00002 0.00000 0.00000 

5d5/25d5/2(0) 5d3/26p3/2(l) 0.00000 0.00000 0.05057 0.07583 0.00000 0.00001 0.00000 0.00000 

5d5/25d5/2(0) 5d5/26p3/2(l) 0.87030 0.75077 -0.15075 -0.26297 0.00100 0.00055 0.87028 0.00003 

6si/26si/2(0) 5d3/26pi/2(l) 0.00000 0.00000 -0.30143 0.03609 0.00006 0.00004 0.00000 0.00000 

6.si/26.si/2(0) 5d:5/26p:!/2(l) 0.00000 0.00000 0.20141 -0.18719 -0.00003 -0.00001 0.00000 0.00000 

6.Si,^().'.i/2(0) 5(/:,/o()/*;j;2(i) I). ()()()()() I). ()()()()() -().."):!."•):! i I). 90911 D.OOOOCi 0.1)1)1)11 0.1)1)1)1)1) 0.1)1)1)1)0 



(a) 

Z'""y 5d3„5d,^{2)-5d,,,6p,^(3) 



\ 

v 



length form 



B \ * \ 

/ V 

Z'^'°'Vo'°'°"°'°W 



70 75 80 85 90 95 100 

Nuclear charge Z 



l(b) 

Z'" \ 5d3„5d,,,(2)-5d,,,6p„,(3) 



^ 0.5 



3.o-o-°-°-°-o^ 



* ^3 



velocity form 
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FIG. 4: Uncoupled matrix elements for 5d3/25dB/2(2)-5d5/26pi/2(3) calculated in velocity and length forms 



7 



TABLE V: Line strengths in length L and velocity V forms 
for even-odd-parity transitions in Re+^. 



RMBPT 



First order 



Level 


Level 


L 


V 


L 


V 


5d'' -"Po 


5d6p 


^1 


0.2632 


0.2312 


0.4952 


0.3749 




5d6p 


1 


0.2424 


0.2277 


0.2673 


0.2027 






^Pi 

1 


0.0645 


0.0651 


0.0532 


0.0391 




5d6p 


^1 


0.0040 


0.0032 


0.0375 


0.0288 


5d^ ^So 




^P 

1 


0.0050 


0.0052 


0.0272 


0.0208 


5d^ ^So 


5d6p 


^P 

1 


0.3261 


0.3437 


0.7024 


0.5236 


5d^ ^Pi 


5d6p 


^Pn 


0.5035 


0.5010 


0.5657 


0.4208 


5d6s ^Di 




3Pn 


1.1657 


1.1792 


1.6949 


1.5356 


5d^ '-"Pi 


5d6p 


^p»1 


0.0080 


0.0081 


0.0242 


0.0179 


5d^ ^Pi 




^'p 

J 1 


0.4143 


0.4127 


0.5065 


0.3774 


5d^ 3 Pi 


5d6p 


^p 

1 


0.0358 


0.0340 


0.0346 


0.0264 


5d6s ^Di 




^'p 

1 


2.1855 


2.2762 


3.0478 


2.7590 








0.0274 


0.0296 


0.0709 


0.0528 


5d^ ''Pi 




^ 2 


0.2803 


0.2573 


0.2959 


0.2251 


5d^ 




^P 

-1 1 


0.1648 


0.1639 


0.1486 


0.1122 




5d6p 




0.1145 


0.1145 


0.1690 


0.1262 


5d6s 




'Po 


1.2140 


1.2595 


1.5450 


1.3925 








0.0380 


0.0361 


0.1735 


0.1284 


5d6s ^Da 






2.4488 


2.7097 


3.3267 


3.0271 




5d6 


1^2 


1.6146 


1.6111 


1.7279 


1.2837 




■5d6 


'Po 

^2 


0.1384 


0.1440 


0.1441 


0.1079 


5d6s ^Ds 




'Fo 

J^2 


0.0020 


0.0027 


0.0026 


0.0024 


5d6s ^Ds 




'Po 
J^2 


2.7197 


2.9674 


3.8172 


3.4676 


5d^ ^^3 


5d6p 




1.6878 


1.5330 


1.9319 


1.4599 


^^3 


5d6p 


^P»3 


0.0528 


0.0563 


0.0308 


0.0219 


5d6s ^Ds 


5d6p 


^P3 


1.7322 


1.8326 


2.3084 


2.0843 


Sd^ ^^3 


5d6p 


^P4 


0.1367 


0.1294 


0.1291 


0.0983 


5d6s ^i33 


5d6p 


^P4 


11.1345 


11.5018 


15.2305 


13.8004 


Sd^ 


5d6p 


^P3 


0.5896 


0.5200 


0.8001 


0.6090 


5d^ 


5d6p 


^^3 


3.6963 


3.6447 


4.0769 


3.0379 


5d2 3^4 


5d6p 


^P3 


0.1442 


0.1561 


0.1513 


0.1095 


Sd^ IG4 


5d6p 


^^3 


0.3118 


0.3027 


0.3497 


0.2616 


5d2 IG4 


5d6p 


^P3 


5.6434 


5.5815 


6.2102 


4.6443 


5d2 ^F4 


5d6p 


^P4 


1.9537 


1.9011 


1.9296 


1.4447 


5d2 IG4 


5d6p 


^P4 


0.1544 


0.1571 


0.1535 


0.1114 



trends of line strengths, oscillator strengths, and transi- 
tion rates. 



A. Transition energies 



In Table Vl , transition energies are compared with re- 
cent measurements ||, |[ ^, |[ 0]. We obtain good 
agreement with the experiment for low- lying 5cP, 5d6s, 
and 5d6p levels. We also we obtain reasonable agree- 
ment (500 -1000 cm~^) for highly-excited 5d5/ levels in 
Re VI, Os VII, and Ir VIII ions. However, we disagree 
substantially (2000 -10000 cm'^) for 5d5/ levels in Pt IX, 
Au X, and Hg XI ions. As mentioned in Refs. jj, ^, D, 
mixing between levels of the 5d5f configuration and the 
core excited configurations with a 5p hole could be very 
important for Pt IX, Au X, and Hg XI ions. Indeed, sub- 



stantial mixing of the 5p^5(P and 5p^5d5f configurations 
was found in Ref. From this, we conclude that the 
present values for 5d5f configurations in Pt IX, Au X, 
and Hg XI are less reliable (~ 1%) than data for Re VI, 
Os VII, and Ir VIII. 



B. Fine structure of the 5d6p triplets 

The fine-structure intervals for the ^P, "^D, and "^F 
terms of 5d6p configuration divided by (Z — 65)^ are 
shown in Fig. ^. The fine structures of these levels do 
not follow the Lande rules even for small Z; the '^P lev- 
els are partially inverted, while the '^D and levels 
show regular ordering of the fine-structure splittings for 
both low and high Z. The unusual splittings are caused 
by changes from LS to jj coupling and by mixing from 
other triplet and singlet states. Comparisons are made 
with experimental data from |, |, |, |, |, § in Fig. |. 
Excellent agreement (0.3% - 3%) is found for the four in- 
tervals 3Pi-3Po, ^D3-^D2, ^^4-3^3, and ^F3-^F2. The 
agreement decreases for the ^P2-^Pi and "^1)2-^-01 inter- 
vals, from 10% for Re+^ to 2% for llg+^°. Experimental 
energies for other ions would be very helpful in confirm- 
ing the Z dependence shown in Fig. |^. 



Line strengths and oscillator strengths in 
Yb-like ions 



Trends of the Z dependence of line strengths and os- 
cillator strengths, are shown in Figs. || and 0. In Fig. ||, 
we illustrate the Z dependence of line strengths for the 
four 5d'^-5d6p transitions calculated (a) by RMBPT and 
(b) in first-order to illustrate once again the impor- 
tance of including second-order matrix elements. Com- 
paring the curves in Fig. ^ and Fig. |^b, we can see 
that the Z dependence is the same for both cases ex- 
cept for the small region Z = 80-88. The singular- 
ities in this region are already present for uncoupled 
dipole matrix elements. The second-order contribution 
for the 5(i3/25d5/2(2)-5(i5/26pi/2(3) uncoupled matrix el- 
ement is shown in Fig. ^ Similar singularities appeared 
for other uncoupled matrix elements between even-parity 
J— 2 states and odd-parity J=3 states. Consequently, 
those singularities also appear in the coupled matrix ele- 
ments and the line strengths shown in Fig. ^. 

The Z dependence is more complicated for transitions 
between even-parity states with J=l and odd-parity 
states with J=0, as seen in Fig. |^, where we give os- 
cillator strengths for four transitions calculated (a) by 
RMBPT and (b) in first-order. The oscillator strength 
curves are not smooth functions of Z, even in first-order. 
The singularities shown in the four curves of Fig. 0b arise 
from the strong mixing between three states of the odd- 
parity complex with J=0. The mixing of 5^5/25^5/2, 
5^3/26^3/2, and 6si/26pi/2 states was discussed earlier in 
Fig. I^a where it was shown that avoided crossings occur 
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TABLE VI: Energy levels (cm ^) in ytterbium isoelectronic sequence. Comparison of RMBPT results with experimental data 
presented in Refs. ||, |, |[ |, |. 



Level 


Re VI 


Os VII 


Ir VIII 


Pt IX 


An X 


HgXI _ 




RMBPT 


n 


RMBPT 


nn 


RMBPT 


Pin 


RMBPT 


n 


RMBPT 




RMBPT 


M 


■'F2 








































8491 


8167 


10604 


10308 


12932 


12672 


15490 


15250 


18284 


18077 


21325 


21175 




14495 


14379 


16117 


15837 


17670 


17254 


19174 


18612 


20638 


19949 


22055 


21255 


bd? ^F4 


15095 


14679 


18409 


18049 


21920 


21615 


25633 


25359 


29550 


29309 


33677 


33481 




16866 


16577 


19852 


19444 


22905 


22436 


26146 


25568 


29568 


28925 


33211 


32513 


Sd^ ^Pi 


19556 


19142 


22624 


22098 


25841 


25222 


29238 


28515 


32832 


32039 


36643 


35802 




26807 


26657 


31496 


31274 


36522 


36253 


41949 


41610 


47818 


47451 


54165 


53918 




28371 


27723 


33932 


33127 


39735 


38878 


45962 


44991 


52588 


51566 


59691 


58835 


5d2 ^So 


51056 


50492 


58851 


57710 


66656 


65022 


74624 


72527 


82848 


80385 


91380 


88897 


5d6s '■^Di 


92656 


92312 


129438 


129244 


169138 


169223 


211579 




256622 




304159 




5d6s ^Da 


94554 


94262 


131488 


131416 


171413 


171485 


214009 




259236 




306918 




5d6s ^D2 


109374 


108709 


148595 


148023 


191018 


190553 


236412 




284633 




335596 




5d6s '-^D-i 


103444 


102712 


142679 


142163 


185041 


184748 


230364 




278518 




329409 






198652 




272504 


282345 


352159 


363633 


437184 




527172 




621600 




5d6p ^F2 


163618 


162134 


209906 


208659 


259044 


258017 


310888 


310064 


365322 


364711 


422250 


422506 


5d6p ^Di 


166538 


166077 


213084 


212862 


262439 


262461 


314456 


314751 


369016 


369624 


426038 


426672 


5dQp ^D2 


176168 


174761 


224886 


223453 


276670 


275872 


331369 


330850 


388865 


388660 


449056 


449335 


5d6p 


177212 


175502 


226212 


224741 


278237 


276994 


333148 


332117 


390829 


390038 


451183 


450786 


5d6p 


182853 


181154 


233874 


232453 


288280 


287122 


345942 


345026 


406770 


406109 


470688 


470408 


5d6p '-^D'i 


186672 


185022 


238186 


236726 


293068 


291792 


351186 


350075 


412435 


411533 


476719 


476134 


5d6p ^Pi 


187037 


185963 


238656 


237682 


293545 


292734 


351634 


350986 


412843 


412413 


477109 


476749 


5d6p ^Po 


189199 


187780 


241043 


239794 


296204 


295142 


354580 


353689 


416091 




480682 


480075 


5d6p ^P4 


195446 


193260 


249334 


247354 


306792 


305005 


367714 


366087 


432023 


430576 


499662 


498310 


5d6p ^P2 


196486 


194539 


250142 


248453 


307362 


305894 


368050 


366803 


432129 


431126 


499533 


498771 


5d6p ^Pa 


197742 


195691 


251254 


249401 


308315 


306697 


368814 


367435 


432644 


431579 


499643 


499109 


5d6p ^Pi 


200506 


200437 


255391 


255246 


313586 


313520 


375052 


375145 


439684 


440120 


507190 


508193 


6s6p ^Po 


261637 




343490 




408884 




460493 


456143 


510225 




558566 




6s6p ^Pi 


265550 




341677 


341844 


394150 


391888 


443224 


438038 


490445 




536100 




6s6p ^P2 


281644 


279156 


333322 


332501 


381777 


379116 


428519 


424403 


473776 




517742 




5d5/ 'G4 


282611 


280575 


333410 


331454 


382478 


380334 


430103 


425349 


476544 




522022 




5d5/ ^'Pa 


282986 


281906 


341072 


340330 


391578 


388229 


440422 


435302 


487930 




534343 




5d5/ ^Ps 


283375 


283559 


334390 


334240 


383305 


380724 


430253 


422093 


475351 




518717 




5d5/ ^Jf4 


284803 


282853 


336142 


334302 


385817 


383508 


434072 


430412 


481126 




527163 




5d5/ =^©3 


286237 


286426 


338249 


337436 


388358 


385282 


436878 


425627 


484128 




530397 




5d5/ =^1/5 


287360 


285346 


339706 


337535 


390378 


387545 


439680 


436230 


487868 




535150 




5d5/ ^^Di 


289100 


290208 


349008 




408218 


406542 


459678 


456053 


509418 




557828 




5d5/ ^D2 


289565 


288853 


348392 


347386 


399983 


397557 


449925 


442048 


498451 




545750 




5d5f ^F4 


292767 


290824 


345648 


343897 


396931 


393669 


446783 


439406 


495294 




542569 




5d5f ^Gi 


294002 


293741 


347287 


346144 


398865 


396001 


449051 


445365 


498219 




546694 




5d5f ^D2 
6s6p ^Pi 


295265 


294688 


354838 


354012 


408411 


405867 


460122 


456798 


510546 




559987 




295457 




357819 


356849 


420585 


421774 


477574 




532630 




586680 




5d5f ^Ds 


295676 


295836 


349683 


349155 


401804 


397003 


452259 


436157 


501264 




549026 




5d5f ^Gs 


295717 


295928 


349829 


349324 


402320 


399128 


453392 


443691 


503255 




552106 




5d5f ^He, 


296143 


292600 


349940 




402394 




453792 




504385 




554394 




5d5f 'P3 


300071 


300794 


356198 


356322 


410808 


408418 


464085 


461735 


516269 




567655 




5d5/ ^P2 


301078 


301178 


371823 




465428 




565480 




671378 




782904 




5d5f ^Pi 


301570 


302600 


362544 


363927 


440737 




534533 




634090 




738793 




5d5/ ^Po 


302331 


303600 


358117 


356587 


434121 




527597 




626622 




730708 




5d5f ^Hs 


303559 


307440 


361850 


364949 


418559 


420291 


473757 


470214 


527612 




580323 




5d5f 'P2 


317422 


317600 


396700 




489065 




589085 




695429 




807602 
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FIG. 5: Fins-structure intervals AE/{Z — 65'^ in cm ^ of bd6p ^P, 5d6p ^D, and 5d6p terms as function of Z. Experimental 
data represented by the * symbol are from [Q, ^, ^ ^ ^, ^, 
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FIG. 6: Line strengths for — 5d6p transitions in Yb-like ions: (a) - RMBPT, (b) - First-order approximation. 
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FIG. 7: Oscillator strengths for transitions from even-parity complex with J=l to odd-parity complex with J=0 in Yb-like 
ions: (a) - RMBPT, (b) - First-order approximation. 
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near Z=76 and 84. 



IV. CONCLUSION 

In summary, a systematic second-order RMBPT study 
of the energies of the 5(P, 5d6s, 6s^, 5d6p, 5d5/, and 6s6p 
states of Yb-like ions has been presented. These calcu- 
lations agree with existing experimental energy data for 
intermediate Z=75 to 80 at the level of 500-1000 cm"^ 
for low-lying 5d6s, and 5d6p levels. They provide a 
smooth theoretical reference database for the identifica- 
tion of lines. 

Also presented is a systematic second-order relativis- 
tic RMBPT study of reduced matrix elements, line 
strengths, oscillator strengths, and transition rates for 
allowed and forbidden electric-dipole transitions in Yb- 
like ions with nuclear charges ranging from Z = 72 to 
100. The dipole matrix elements include retardation and 
correlation corrections from Coulomb and Breit interac- 
tions. Both length and velocity forms of the matrix el- 
ements are evaluated, and 5 % differences, caused by 
the non-locality of the starting DF potential, are found 
between the two forms. 

We believe that our results will be useful in analyz- 
ing existing experimental data and planning new exper- 
iments. There remains a paucity of experimental data 
for many of the higher ionized members of this sequence, 
both for term energies and for transition probabilities and 
lifetimes. Additionally, matrix elements from the present 
calculations will provide basic theoretical input for calcu- 
lations of reduced matrix elements, oscillator strengths, 
and transition rates in three-electron Lu-like ions. 
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APPENDIX A: THE SECOND-ORDER 
DOUBLE-EXCITATION CONTRIBUTION 



Here, we consider the second-order double-excitation 
contribution to explain singularities in the Z dependence 
of energy levels given in Figs. |2| and |^. A typical contri- 
bution from one of the double-excitation diagram for the 



second-order interaction energy has the form 
Ai {vw; v'w') oc 2_] X/ — , ^ \ : , (Al) 



mn kk' 



where Xk' (vwmn) are products of angular coupling coef- 
ficients and Slater integrals pT| |, e,„ is the single-particle 
energy for state to, and v designates a single-particle va- 
lence state (n^lvjv)- For the case of a [Xe]4/^'* core, 
sums over m and n include 5^3/2, 5^5/2, 5/5/2, 5/7/2, 
5(77/2, 5(79/2, and all states with principal quantum num- 
ber n >5. We must, however, exclude terms with pairs 
(rnn) that are included in model space. Therefore, we re- 
move the pairs {5d5d), (5d6s), and (6s6s) since they are 
in the even-parity model space and pairs (5c?6p), {5d5f), 
and (6s6p) in the odd-parity model space. 

As an example, let us consider the singularity near 
Z—85 in the curve describing the energy of the 6s^ 
level in Fig. The denominator in Eq. (Al) 

is D ~ 2ees ~ fm ~ for the matrix element 
^;(6si/26si/2; 6S1/26S1/2). For an ion with Z=85, D be- 
comes very small when (TOn)=(5(i3/26(i3/2) since e6si/2~" 
9.488404 a.u., e5d3/2=-12. 069226 a.u., and £6^3/2=- 
6.932272 a.u., giving i:»=0.02469. This one term dom- 
inates the matrix element ^;(6si/26si/2; 6S1/26S1/2) and 
increases its size by a factor of six in comparison with 
values for neighboring Z, explaining the singulatity at 
Z=85 in Fig. gb. The explanation of the sharp features 
in the curves describing energies of the 6s]^/26pi/2 and 
5^5/25/5/2 levels in Fig. ^ is similar. 



APPENDIX B: THE SECOND-ORDER DIPOLE 
MATRIX ELEMENT 

A typical contribution from one of the second-order 
correlation corrections to the dipole matrix element 
{vw{J) — v'w' {J') has the form 

- v'w'ij')] oc y: , 

6j -f- €nj Cy' ^w' 

I 

where Zi^, is a single-electron dipole matrix element. In 
the sum over i, only terms with vanishing denominators 
are excluded. For the 5(i5/25d3/2(2) — 5(i5/26pi/2(3) tran- 
sition, we obtain 

Z(™") [5^5/25^3/2(2) - 5^5/26^1/2(1)] 

^(»,5rf5/2)Afc(5d5/26pi/2 5d3/2z) 



+ e(5c?3/2) - £(5^5/2) - e(6pi/2) ' 



For the case 1=5/7/2 and Z~8A, the denominator be- 
comes 

e^ + e(5d3/2) - e(5d5/2) - e(6pi/2) = 0.003. 

This one term dominates the entire matrix element. For 
the 5(^5/25(^5/2(0) — 5(i3/26pi/2(l) transition, the denom- 
inator becomes very small at Z=88 when 1=5/7/2. 

e^ + e(5(i5/2) - £(5^3/2) - e(6pi/2) = -0.056. 
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Again, this single term dominates the matrix clement. 

A typical contributions from the second-order RPA 
correction for dipole matrix element {vw{J) — v'w'{J') 
has the form 



Zr''^>[vw{J) - v'w'iJ')] o^Yl 



+ 



where the index b designates core states and n des- 
ignates an excited state. For the special case of the 
5rf5/25rf3/2[2] — 6pi/25rf5/2[l] transition, we obtain 



oc 



EE 

n b 



Z{b,n)Xk{5d3/2b 6pi/2n) 
e„ + e(5d3/2) - e(6pi/2) - £6 



In the case of b=5p^/2 and n=5ds/2 for nuclear charge 
Z=88, the denominator becomes 

£(5^3/2) + e(5d3/2) - e(6pi/2) - e(5p3/2) = -0.278. 

As before, the small value of the denominator leads to 
an anomalous increase in the size of the RPA matrix 
element. 



ZpP^)[5d3/25ci5/2(2) - 6pi/25d5/2(3)] 
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